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Abstract 

An implementation of the strong-permittivity-fluctuation theory (SPFT) is presented in order to 
estimate the constitutive parameters of a homogenized composite material (HCM) which is both cu- 
bically nonlinear and anisotropic. Unlike conventional approaches to homogenization, the particles 
which comprise the component material phases are herein assumed to be small but not vanish- 
ingly small. The influence of particle size on the estimates of the HCM constitutive parameters 
is illustrated by means of a representative numerical example. It is observed that, by taking the 
nonzero particle size into consideration, attenuation is predicted and nonlinearity enhancement is 
somewhat diminished. In these respects, the effect of particle size is similar to that of correlation 
length within the bilocally-approximated SPFT. 
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1 INTRODUCTION 

Formalisms based on the strong-permittivity-fluctuation theory (SPFT) have been developed to es- 
timate the constitutive parameters of homogenized composite materials (HCMs), within the realms 
of both linear [1, 2, 3, 4, 5, 6] and weakly nonlinear [7, 8, 9, 10] materials. The SPFT approach 
to homogenization has the advantage over more conventional approaches, such as those named 
after Maxwell Garnett and Bruggeman, that the distributional statistics of the component material 
phases is generally better taken into account [11]. For example, the SPFT is most commonly imple- 
mented at the level of the bilocal approximation, wherein the distributional statistics are described 
in terms of a two-point covariance function and its associated correlation length. Thereby, the 
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bilocally-approximated SPFT predicts attenuation in the HCM, even when the component phase 
materials are nondissipative. 

In the SPFT approach to homogenization, and likewise the Maxwell Garnett and Bruggeman ap- 
proaches, the electromagnetic responses of the particles which make up the component material 
phases are represented by electrically-small depolarization regions [11]. Often, these depolariza- 
tion regions are taken to be vanishingly small [12, 13]. The spatial extent of the component phase 
particles is thereby neglected. However, the importance of the nonzero spatial extent of the com- 
ponent phase particles has been underlined by several studies pertaining to linear isotropic HCMs 
[14, 15, 16, 17, 18, 19]. Furthermore, versions of the SPFT which accommodate depolarization re- 
gions of nonzero volume were recently developed for linear anisotropic [20] and bianisotropic 3 [21] 
HCMs. While these studies take into consideration the nonzero spatial extent of the component 
phase particles, it is essential that the particle sizes are much smaller than the electromagnetic 
wavelengths, in order for the assembly of component material phases to be regarded as an effec- 
tively homogeneous material [11]. Accordingly, at an optical wavelength of 600 nm, for example, 
component phase particles with linear dimensions less than approximately 60 nm are envisaged. 

In the present study, we extend the consideration of depolarization regions of nonzero volume 
into the weakly nonlinear regime for anisotropic HCMs. Our starting point is the SPFT for cu- 
bically nonlinear HCMs which incorporates vanishingly small depolarization regions [10]. In §2 
the component material phases and their statistical distributions are described. The linear and 
weakly nonlinear contributions to the depolarization dyadics, arising from depolarization regions 
of nonzero volume, are presented in §3. The corresponding bilocally-approximated SPFT estimate 
of the HCM permittivity dyadic is given in §4; and a representative numerical example is used to 
illustrate these results in §5. Finally, a few closing remarks are provided in §6. 

In the notation adopted, single underlining denotes a 3 vector whereas double underlining denotes 
a 3x3 dyadic. The inverse, determinant, trace and adjoint of a dyadic M are represented as M^ 1 , 
det ( M ) , tr (M) and M ad i , respectively. The permittivity and permeability of free space are 
written as e and fi . 

2 HOMOGENIZATION PRELIMINARIES 

The homogenization of two component material phases, namely phase a and phase b, is considered. 
Each component phase is an isotropic dielectric material; and, in general, each is cubically nonlinear. 
Thus, the permittivities of the component material phases are expressed as 

Q = Qo + XI \Ke ! 2 > (Z = a,b), (1) 

with e£Q being the linear permittivity, xe the nonlinear susceptibility, and E_ g is the electric field 
developed inside a region of phase £ by illumination of the composite material. The assumption 
of weak nonlinearity ensures that j qo | 3> | Xi I \Ke \ 2 ■ Notice that nonlinear permittivities of the 
form (1) describe electrostrictive materials which can induce stimulated Brillouin scattering [22]. 

The component material phases a and b are made up of ellipsoidal particles. The particles all have 
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the same shape and orientation, as specified by the shape dyadic 



U= diag{U x ,U y ,U z ), (U x ,U y ,U z >0), (2) 

— yu x u y u z 

which parameterizes the particle surface as 

r e (6,^=r ] g>r(e,4>), (3) 

where f (9, <p) is the radial unit vector specified by the spherical polar coordinates 6 and 4>. The size 
parameter rj provides a measure of the linear dimensions of the ellipsoidal particles. It is assumed 
that r] is much smaller than the electromagnetic wavelengths, but not vanishingly small. 

The component phase particles are randomly distributed throughout a region of volume V, which 
is partitioned into the disjoint regions of volume V a and V& containing phase a and 6, respectively. 
Thus, the component phase distributions are characterized in terms of statistical moments of the 
characteristic functions 

(I, r£ V e , 

Ml) = { (£ = a,b). (4) 

,0, r $ V e , 

The first statistical moment of <3?^ delivers the volume fraction of phase £, i.e., ( $^(r) ) = fy. Plainly, 
fa + fb = 1- The two-point covariance function which constitutes the second statistical moment of 
is taken as the physically-motivated form [23] 

( (Ml))(Ml')), \UT l • (r-r')l >L, 

{Ml)Ml')) = { (5) 

I (ML)), \UT l ' (L-r')\<L, 

where L > is the correlation length. Within the SPFT, the estimates of HCM constitutive 
parameters are largely insensitive to the specific form of the covariance function, as has been shown 
by comparative studies [8, 24]. 



3 DEPOLARIZATION DYADIC 



Let us focus our attention on a single component phase particle of volume V e , characterized by 
the shape dyadic U and size parameter r/. Suppose that this particle is embedded in a comparison 
medium. In consonance with the ellipsoidal geometry of the component phase particles and the 
weakly nonlinear permittivities of the component material phases, the comparison medium is a 
weakly nonlinear, anisotropic, dielectric medium characterized by the permittivity dyadic 

Icm = IcmO + ^ cm I —HCM ? = dia g (4o> e L ' e cm ) + dia g (Xcm, Xlm, Xcm) \ElhCM ? \ (6) 

where Elhcm denotes the spatially-averaged electric field in the HCM. The eigenvectors of § cm are 
aligned with those of U_. 
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The depolarization dyadic [12] 



DM = f G (r) d 3 r 



(7) 



provides the electromagnetic response of the ellipsoidal particle embedded in the comparison 
medium. Here, the dyadic Green function of the comparison medium, namely G (r), satisfies 
the nonhomogeneous vector Helmholtz equation 

(VxVxI-wVolJ -a cm (r-r!)=iu^S(r-r!)L. (8) 
An explicit representation of G_ cm {r) is not generally available [25], but its Fourier transform, 

icJD = I a c Jr-) exp(-ig • L) d 3 r, (9) 

Jr_ 

may be deduced from (8) as 

£ante) = - i(J Vo («x gx| + o; 2 /xo| cm ) 1 . (10) 
By combining (7), (9) and (10), 



D 



_n_ f ]_ 

2vr2 J g 



sin (qrj) 



qrj 



cos (qrj) 



G(U- 1 'q)d 3 q 



is obtained, after some simplification [12, 26]. 



(11) 



3.1 Depolarization contributions from regions of nonzero volume 

As in [20], we express the depolarization dyadic as the sum 

D z = Dj F=o + Dji >0 , 

where the dyadic 



2vr2 / q* 



sin (qrj) 
qrj 



— cos(qrj) 



lim GJSL^'q) 

q— >oo — cm — — 



d 3 q 



(12) 



(13) 



represents the depolarization contribution arising from the region of vanishingly small volume 

lie 

sin (qrj) 



limV" e , whereas the dyadic 

— >0 



n v>o = _n_ 



lq 2 



qrj 



— cos(qrj) 



lim G OJ- 1 - q) 

q^oo =™ - - 



<Fq (14) 



provides the depolarization contribution arising from the region of nonzero volume I V e — \\m.V e 

ij—fO 
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Depolarization dyadics associated with vanishingly small regions have been studied extensively 
[13, 26]. The volume integral (13) reduces to the ^-independent surface integral [12] 



2tt 



with 



A = diag 



o Je=o 



sin 2 9 cos 2 . 



tr ( e • A 

--cm = 



A 



sin dOdcfi, 



sin 2 9 sin 2 6 cos 2 > 



(15) 



(16) 



An elliptic function representation for iy? =0 is available [27] (which simplifies to a hyperbolic 
function representation in the case of a spheroidal depolarization region [12]), but for our present 
purposes the integral representation (15) is more convenient. 

Depolarization dyadics associated with small regions of nonzero volume have lately come under 
scrutiny for anisotropic [20] and bianisotropic [21] HCMs. As described elsewhere [10, 20], by the 
calculus of residues (14) reduces to 

DV >o = _±_w(r}), (17) 



where the dyadic function 

tttt / \ 3 f 2w r sin9 

WW = v s 

J 6=0 Je=o 



3A 



3 (k + — K. 



27? 



+ i Kl — KZ. 



g_ + iu iJ, — K-) §_[ d9d<fi 



is introduced. Herein, the dyadics 



a = 



v -, , v tr(e°# -A) - \ti(e ad i) ti ( A) 

26 -trfe )l -A-trfe • a) I - V=cm =' . L ^ ' K = J 

= cm \ = cm I = = \ = cm = I = I . \ 

tr e -i 

y = cm = J 



P = i 



adj 
cm 



det 



\ = cm J 



trfe • A) 

\ = cm = J 



A 



and the scalars 



A = ^Jt 2 B - At A tc, 
2 -t B ± A 



K± = fi UJ 



2t c 



with 



tA = det 
t B = tr ( 



( e ) 

\ = cm J 
£ adj . 



tr tr (A) 



t c = tr(e =cm - A) tr (A) 



(18) 

(19) 
(20) 

(21) 
(22) 

(23) 



Often the approximation D_ ss jyn =0 is implemented in homogenization studies [11]. However, 
studies of isotropic [14, 15, 16, 17, 18, 19], anisotropic [20] and bianisotropic [21] HCMs have 
emphasized the importance of the nonzero spatial extent of depolarization regions. 
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3.2 Linear and weakly nonlinear depolarization contributions 



We exploit the fact that the comparison medium permittivity (6) is the sum of a linear part and a 
weakly nonlinear part to similarly express 



n z = n z0 +n zl \E HCM \ 2 



RT° + =T° + (=r° + =T°) I - 



HCM I > 



where 



^=Rf°+Rf°\E HCM f. 



(24) 



(25) 



The linear and weakly nonlinear contributions to D_ v have been derived earlier [10]; these are 



»7=0 



— Airitu 



2n rn 



J 6=0 



■A 



Q = _ 1 

4vriu; J(j>=0 



f'TV 

Je=o 



tr (i cm0 -4) 

tr( X -A) 

\ = cm — / a 



sin 9 d9d(j), 



tr 



(icmo'i) 



2A 



> sin d6 d<t>. 



(26) 



(27) 



The linear and weakly nonlinear contributions to D v>0 — and, equivalently, Wjrj) — follow from 
corresponding contributions for an expression analogous to (18) which crops up in the bilocally- 
approximated SPFT [5, 10]. Thus, we have 



m(v)=m (v)+mi(v)\E HCM \ 



with 



and 



K (v)=v 3 



2tt r-w 



sm I 



J 9=0 3 A 



de d4> 



(28) 
(29) 



(30) 



where 



^(l) = 3 ^ 0+ 2r? K ° - + * ( 'Ml- - '-()- 



1 1 



(31) 
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The dyadics a Q and /3 , and scalars kq± and Ao, herein represent the linear parts of their counter- 
part dyadics a and 0, and scalars k± and A, as per [10] 



^=cm0 ^( = cm()) = = ^ (=cm0 =) = 



tr ( e adj n • A 

:cm0 = 



tr 



(lcmo'4) 



= e adj 

= o = cmO 



det 



( = cmo) 



«0± = Mo^ 



te (l cm 0-i) 

2 —ten ± A 



A, 



Ao = yj t 2 m - 4:tAotco, 



tr ri tr (4) 



with 

t A0 = det(| cm0 ) 

ico = tr( |cm0 .4)tr(4) 
Moreover, the weakly nonlinear contributions to a, /?, k± and A are provided as [10] 

tBitco — tBotci 



Ai = 



2x 



tcotr(g cm0 -A) 



tr 



(x )/ 

\ = cm / — 



A-tr( X 'A)i, 

— \ = cm — / — 



^ _ tBitco — tsotci 

tcotr[ gem0 -A)- 

_ ( -tm ± Ai ) - 2fgi k 0± 
2tco 

tBptBl - 2 (tAl^CO + tAptd ) 
An 



with 



and 



, x y z i x y z i x y z 

tAl = Xcm e cm0 e cm0 + e cm0 Xcm £ C m0 + e cm0 e cm0 Xcm 

t B1 =tv(T'A) - [tr(X) tr(^)] 
*Ci =tr(A) tr fx 'A) 

— \ = cm — / 

T = (\\s\p ( v v f z + f y v 2 v 2 f 1 -I- <= 2 v x v x f y A-e x v y \ 
— u a 6 V Acm c cm0 ^ c cm0 Acmi Acm c cm0 ^ c cm0 Acm) Acm c cm0 ^ cmO Acm / ■ 



4' 

(32) 
(33) 

(34) 
(35) 



(36) 



(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



4 SPFT ESTIMATE OF HCM PERMITTIVITY 

Now that the linear and nonlinear contributions to the depolarization dyadic have been established 
for depolarization regions of nonzero volume, we can amalgamate these expressions with the SPFT 
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for weakly nonlinear anisotropic dielectric HCMs — which is presented elsewhere [10] — and thereby 
estimate the HCM permittivity. 

As a precursor, an estimate of permittivity dyadic of the comparison medium must first be com- 
puted. The Bruggeman homogenization formalism (which is, in fact, equivalent to the lowest-order 
SPFT [5]) is used for this purpose. Thus, e is found by solving the nonlinear equations 



faZ ai +f b 2L br g, 



(j = 0,1), 



(43) 



where 



^i = ~ iuj ( e ^z-i cm0 ) 



(£ = a,b), (44) 



are the linear and nonlinear parts, respectively, of the corresponding polarizability dyadics. Herein, 
1 



det 



(e») 



diag ^r^rfg + T y iQ v z n , r^rfo + Ti Q T x n , r^rf + r^ rf 1 ^ - p£]_ t l 



with 



„ _ -px -pV -pz , -px -pV -pz . -px -pV -pz 
Pi — L £0 L eo L £1 "r i ^0 i £l i £0"r i £l i ^0 i £0' 

are expressed in terms of components of the dyadics 



(45) 



(46) 



= £0 = ' =0 



I + iuD n • ( e eo l-g cm0 



dia g (r? ,r^ , r| ) 
fi ' (gexiL-z cm ) + £i * (^oI-| cm0 ) 

and the local field factor is estimated by [28] 



diag^r^rfj 



(47) 



tr 



(48) 



Estimates of the i cm0 and x may be straightforwardly extracted from (43) by recursive schemes; 
see [10] for details. 

Finally, the bilocally-approximated SPFT estimate of the HCM permittivity dyadic, namely 

In = Ino + K n I ^hcm I 2 = diag (4 , e y no , e z no ) + diag ( x fi, x% Xn) I ^tfCM | 2 , (49) 
is given as [5] 



l™ = e 



= i l = cm 



- —Q' 1 • s 

iw= — u 

-— (V 1 -s l + n-s 

\= — 1 — — u 



(50) 
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Herein, the linear and nonlinear parts of the mass operator are represented, respectively, by the 
dyadics 



S„ = 



— Aniuj 

fafb 



> ; 



and the dyadic 



n 



det ( Q q 



with 



v — QqQqQi + Q0Q1Q0 + QiQqQo, 
is expressed in terms of the components of 



£ =I + Io -£ = dia g^o> Qo> Qo) 



(51) 



diag ( QXQl + QlQl, QlQl + QgQf, Qjf Qg + Q^Qf ) - i/g" 1 ] , (52) 



(53) 



(54) 



5 NUMERICAL STUDIES 



The SPFT estimates (50) of the HCM linear permittivity and nonlinear susceptibility are rep- 
resented by mathematically complicated expressions. In order to discern the influence of the 
size parameter rj, parametric numerical studies are called for. To this end, we investigate the 
following representative example of a homogenization scenario. Let component phase a be a cu- 
bically nonlinear material with linear permittivity e a Q = 2e and nonlinear susceptibility Xa = 
9.07571 x 10~ u e„ m 2 \- 2 (= 6.5 x 10~ 4 esu); and component phase b be a linear material with 
permittivity e& = e&o = 12e . The eccentricities of the ellipsoidal component phase particles are 
specified by U x = 1, U y = 3 and U z = 15. These choices of parameter values facilitate direct 
comparisons with a previous investigation in which the effects of the size parameter rj were not 
included [10]. Results are presented for an angular frequency of to = it x 10 15 rad s _1 (equivalent to 
a free-space wavelength of 600 nm). 

We begin with the relatively straightforward case where neither the size parameter nor the correla- 
tion length is taken into account; i.e. rj = L = 0. In this case, the SPFT estimates of the constitutive 
parameters are equivalent to those of the conventional Bruggeman formalism for weakly nonlinear, 
anisotropic, dielectric HCMs [28]. In Fig. 1, the HCM linear and nonlinear constitutive parameters 
are plotted against volume fraction f a . The HCM linear permittivity parameters £qq' z uniformly 
decrease from e^o at f a = to e a o at f a = 1. In contrast, the HCM nonlinear susceptibility pa- 
rameter xfii an d to a lesser extent xfii exceeds the nonlinear susceptibility of component phase a 
for a wide range of values of f a . This nonlinearity enhancement phenomenon, and its potential 
for technological exploitation, have been reported on previously for both isotropic [7, 9, 29, 30] 
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Figure 1: The HCM relative linear permittivity and nonlinear susceptibility parameters plotted 
against f a , calculated for 77 = L = 0. Key: eQ /e and xfi/Xa dashed curves; e^ / e o an ^ Xn/Xa 
broken dashed curves; and £Q / e o ano - Xn/Xa solid curves. Component phase parameter values: 
e a0 = 2e , X a = 9.07571 x 1Q- I2 e m 2 \- 2 , e b = e m = 12e„, U x = 1, U y = 3 and C/ 2 = 15. 




Figure 2: Real and imaginary parts of the HCM linear permittivity and nonlinear susceptibility 
parameters plotted against r] (in nm), calculated for L = and f a = 0.3. Key: e^ 0r and x^r 
dashed curves; e^ 0r and Xs> broken dashed curves; and e^ 0r and Xn r solid curves. Component 
phase parameter values as in Fig. 1. 

and anisotropic [10, 28] HCMs. The anisotropy reflected by the constitutive parameters, and the 
nonlinearity enhancement, stems from the ellipsoidal geometry of the component phase particles. 

How does the size parameter r] influence the estimates of the HCM constitutive parameters? To 
answer this question, we fix the volume fraction at f a = 0.3 and calculate the HCM constitutive 
parameters for < rj < 20 nm with L = 0. The presentation of results is aided by the introduction 
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of the relative constitutive parameters 



e no \ e no\ v =L=o 



~Q0r 



Xn [Xn 1^=^=0 

Xa 



(n = x,y,z), 



(55) 



which measure the difference between the SPFT estimates calculated for ij,L/0 and r] = L = 0. 
The results are plotted in Fig. 2. It is notable that the HCM constitutive parameters have nonzero 
imaginary parts whereas the component material phases are specified by real-valued constitutive 
parameters. As previously described for linear HCMs [20], the presence of nonzero imaginary parts 
for the HCM constitutive parameters may be attributed to radiative scattering losses associated 
with the nonzero size of the component phase particles. Plainly, increasing the size parameter r] has 
the effect of increasing the real and imaginary parts of the HCM linear permittivity, but decreasing 
the real and imaginary parts of the HCM nonlinear susceptibility. In fact, the influence of the size 
parameter is very similar to the influence of the correlation length, as has been noted for linear 
HCMs [10]. 



Im e QOr 




Re xtt 




L (nm) 




L (nm) 




L (nm) 



Figure 3: Real and imaginary parts of the HCM linear permittivity and nonlinear susceptibility 
parameters e^ 0r and xfi r plotted against L (in nm) and rj/L, calculated for f a = 0.3. Component 
phase parameter values as in Fig. 1. 
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Fig. 2 reveals that by taking into consideration the nonzero size of the component phase particles — 
but not the correlation length — the predicted nonlinear ity enhancement is somewhat diminished. 
We now consider the estimates of the HCM constitutive parameters when both the size param- 
eter and the correlation length are taken into account. In Fig. 3, the HCM relative constitutive 
parameters are plotted against both L and rj/L with the volume fraction fixed at f a = 0.3. Only 
the results for e^Q and xfi are presented; the corresponding plots for 6qq and Xq Z are similar. It 
may be observed in Fig. 3 that the effects of i] and L are cumulative insofar as the increase in the 
real and imaginary parts of e^'q' 2 , and the decrease in the real and imaginary parts of Xq V ' Z \ which 
occur as rj increases, become steadily more exaggerated as L increases. 

6 CONCLUDING REMARKS 

The size of the component phase particles can have a significant bearing upon the estimated consti- 
tutive parameters of weakly nonlinear anisotropic HCMs, within the bilocally-approximated SPFT. 
Most obviously, by taking nonzero particle size into consideration, attenuation is predicted and the 
degree of nonlinearity enhancement is somewhat diminished. In respect of both of these effects, the 
influence of particle size is similar to the influence of correlation length. Furthermore, the effects of 
particle size and correlation length on both the linear and nonlinear HCM constitutive parameters 
are found to be cumulative. 
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